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ABSTRACT 

It is a well-established empirical fact that the surface density of the star formation rate, Esfr, 
strongly correlates with the surface density of molecular hydrogen, Eh 2 , a t least when averaged over 
large (~ kpc) scales. Much less is known, however, if (and how) the Esfr — Eh 2 relation depends on 
environmental parameters, such as the metallicity or the UV radiation field in the interstellar medium 
(ISM). Furthermore, observations indicate that the scatter in the Esfr - Eh 2 relation increases rapidly 
with decreasing averaging scale. How the scale-dependent scatter is generated and how one recovers a 
tight ~ kpc scale Esfr — Eh 2 relation in the first place is still largely debated. Here, these questions 
are explored with hydrodynamical simulations that follow the formation and destruction of H2 , include 
radiative transfer of UV radiation, and resolve the ISM on ~ 60 pc scales. We find that within the 
considered range of H2 surface densities (10-100 Mq pc -2 ), the Esfr — Eh 2 relation is steeper in 
environments of low metallicity and/or high radiation fields (compared to the Galaxy), that the star 
formation rate (SFR) at a given H2 surface density is larger, and the scatter is increased. Deviations 
from a "universal" Esfr — Eh 2 relation should be particularly relevant for high-redshift galaxies or 
for low-metallicity dwarfs at z ~ 0. We also find that the use of time-averaged SFRs produces a large, 
scale-dependent scatter in the Esfr — Eh 2 relation. Given the plethora of observational data expected 
from upcoming surveys such as ALMA, the scale-scatter relation may indeed become a valuable tool 
for determining the physical mechanisms connecting star formation and H2 formation. 
Subject headings: galaxies: evolution — methods: numerical — stars: formation 



1. INTRODUCTION 



In a seminal paper, ISchmidtl (|1959f ) constructed a 
closed-box model of gas consumption and star forma- 
tion that relies on the basic assumption of a polyno- 
mial relationship between (total) gas surface density E gas 
and star formation rate surface density Esfr- This 
model was able to satisfy simultaneously a number of 
observational constraints, such as the initial luminos- 
ity function of main-sequence stars, the luminosity func- 
tion of white dwarfs, or the relatively constant surface 
density of atomic hydrogen (HI). While the first stud- 
ies focusse d on the relation between neutra l hydrogen 
and SFR (jSanduleakl Il969t iHartwickl Il97lh . the com- 
bination of measurements of Ha, HI and CO emission 
lines allowed for a direct test of the Schmidt relation, 
£-Sfr — E gas , and a precise measurement of its expo- 
nent ([Kennicu ttl fl989L [19980 ). Initially, it was assumed 
that E gas would determine Esfr (e.g., via gravitational 
collapse). However, measurements of azimuthally av- 
eraged gas and SFR profiles showed that SFRs cor- 
relate better with the molecular hydro gen (H2) com- 
pone n t than with the t otal gas density ()Wong fc Blitzl 
l2002t IBigiel et all I2008T ). In fact, recent observational 
and theoretical works demonstrate that the steepening 
of the Esfr — E gas relation at low gas surface densi- 
ties coincides with the transition of atomi c to molecu- 
lar hydrogen (jRobertson fc Kravtsov! 120081 IBigiel et~aT1 
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[20081 iKrumholz et al.ll2009bl iGnedin fc Kravtsovll20H . 
The shape of the Esfr — E gas relation is also predicted 
to evolve strongly with redshift due to the build-up of 
metallicity in the interstellar medium (ISM) over cos- 
mic history and the importance of dust in the forma- 
tion of H2 and its shielding from Lyman- Werner radia - 
tion (jKrumholz et al.ll2009at IGnedin fc Kravtsovll2010t ). 
In contrast, the Esfr — Eh 2 relation is often assumed 
to evolve little and be relatively insensitive to changes 
in metallicity and interstellar radiation field, although 
this has not yet been confirmed observationally. The as- 
sumption on which this "universality" is based is that 
the efficiency with which clouds of molecular hydrogen 
convert their H2 into stars is not a strong function of 
the average ISM metallicity or the interstellar radiation 
field, at least under con ditions typical for spiral galaxies 
(jKrumholz fc T an 2007). The scale at which this conver- 
sion takes place is the scale of (giant) molecular clouds, 
i.e., 100 pc or less. However, there are a couple of com- 
plications. First, there is growing observational evidence 
suggesting that the scatter in the Esfr — Eh 2 relation in- 
creases if one goes to smaller and smaller scales (see, e.g., 
lOnodera et"aJ]|2010t ISchruba et al.ll2010D . Taken at face 
value, this seems to contradict a tight small-scale cou- 
pling between molecular hydrogen surface density and 
star formation. Second, the Esfr — Eh 2 relation is typ- 
ically measured on ~ kpc scales and the spatial aver- 
aging may lead to changes in the slope, inte rcept and 
scatt er compared with those on small scales ([Kravtsovl 
2003). Third, observationally determined SFRs are time- 
averaged over the effective lifetime of the specific star 
formation tracer and may thus differ from instantaneous 
SFRs. 

A straightforward observational check of the univer- 
sality of the Esfr — Eh 2 relation is difficult, first and 
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TABLE 1 

Overview of the simulations discussed in this work. 



n c (cm 3 ) Z/Zq Umw Resolution (pc) No. 



50 0.1, 0.3, 1 0.1, 1, 10, 100 65 12 

10 3 , 10 6 1 0.1 65 2 

10 3 , 10 6 0.1 100 65 2 

50 1 1 32 1 

50 1 1 125 1 



stars at a given H2 density independent of the scatter 
width a. This implies, however, that star formation ef- 
ficiencies derived from the mean (or median) relation in 
the log p* — logpn 2 plane, i.e., for X = 0, differ from the 
parameter csfr by the factor e~2 CT . This effect should 
be considered when, e.g., star formation efficiencies of 
scatter-free numerical simulations are calibrated against 
observations on sub-kpc scales. 

On which grounds would we actually expect to see 
any dependence of slope, intercept and scatter of the 
EgpR— Sjj 2 relation on environmental parameters such as 
metallicity or interstellar radiation field? The H 2 abun- 
dance is strongly affected by the amount of dust shielding 
from the UV radiation and, consequently, a lower metal- 
lically and/or larger radiation field will increase the re- 
quired density for H2 (and consequently stars) to form. 
We will show that a non-linear relation between nn 2 and 
the SFR on small scales can have a significant impact on 
the slope, intercept, and scatter of the Ssfr — £r- 2 rela- 
tion measured on large (~ kpc) scales. Another impor- 
tant, and so far often neglected quantity, is the scatter in 
the Ssfr — Sh 2 relation. While some scatter may be due 
to observational measurement uncertainties, it is clear 
that any environmental dependence of the Esfr ~ Sjj 2 
relation will translate into a galaxy-to-galaxy variation 
and, in combined data sets, to scatter. Furthermore, the 
observed Esfr — Sh 2 relation is measured on large scales 
(spatial averaging) using time-averaged SFRs. The aver- 
aging may induce a scale-dependence of the scatter. 

The layout of the paper is as follows. In fj^lwe briefly 
describe the setup of our numerical experiments. We 
then show in §3.11 the predicted dependence of the slope, 
intercept, and scatter of the Esfr — Xh 2 relation on 
metallicity and interstellar radiation field. The scale de- 
pendence of the scatter and the propagation of intrinsic 
scatter from 100 pc to kpc scales is studied in £13.21 We 
discuss our findings in £0] and conclude in <J5] 

2. SIMULATIONS 

A detailed descriptio n of the set of performed simu- 
lations can be found in iGnedin fe Kravtsovl (|2011l ). All 
simulations are run with the Eulerian hydrodyn amics + 
A-body code ART (jKravtsov et al.lll997l 12002ft . which 
uses an adaptive mesh refinement (AMR) technique to 
achieve high spatial resolution in the regions of inter- 
est (here: regions of high baryonic density). First, we 
ran an initial cosmological, hydrodynamical simulation 
down to z = 4. This simulation follows a Lagrangian re- 
gion that encloses five virial radii of a typical galaxy 
(halo mass ~ 10 12 M Q at z = 0) within a box of 6 co- 
moving Mpc//i. The mass of dark-matter particles in the 
high-resolution Lagrangian patch is 1.3 x 10 6 Mq and the 
spatial resolution is 65 pc at z = 3 in physical coordi- 
nates. We adopt the following cosmological parameters: 
^matter = 0.3, Cl A = 0.7, h = 0.7, f^baryon = 0.043, and 
as = 0.9. This initial, fully self-consistent simulation 
is consequently continued for additional ~ 600 Myr be- 
fore it is analyzed, but now with metallicities and UV 
fields fixed to a specific, spatially uniform value. At this 
time, the mass of the simulated halo is s» 4.2 x 10 11 M@. 
We have run a grid of simulations with three different 
metallicities Z = 0.1,0.3,1.0 (in units of Z Q = 0.02) 
and four different values of the interstellar radiation field 
Umw — 0.1,1,10,100. The parameter Umw = J/Jmw 



Note. — The columns denote: (1) the minimum density of 
molecular clouds that form stars n c , see eqution JSJ, (2) the metal- 
licity, (3) the radiation field, (4) the spatial resolution, and (5) the 
number of the simulations discussed in this work. 

foremost because the direct detection of molecular hy- 
drogen is challenging. Tracer molecules such as CO or 
HCN are ty pically used in stead to infer the H2 column 
density (e. g. . I0montll2007l ) . Mapping the line intensity of 
tracer molecules to the H2 column density is obfuscated 
by the fact that the dependence of the conversion factor 
on ISM properties, e.g., metallicity or interstellar radi- 
ation field, is not well understood (|Glover fc Mac Lowl 
1201 ll) . In addition, rad iative transfer effects need to be 
carefully modeled (e.g., iNaravanan et al.ll2008l 120111 ). 

Numerical simulations offer a different route to study- 
ing the Ssfr — Sh 2 relation. What are their require- 
ments? First, the numerical code needs to follow self- 
consistently the formation and destruction of H2. This 
implies a resolution of 100 pc or better, an implementa- 
tion of cooling down to a few tens of Kelvin and, also, 
radiative transfer of the Lyman- Werner bands (at least 
in some approximate form), in order to correctly capture 
the impact of the interstellar radiation field on the H2 
dissociation. 

Secondly, the code needs a recipe for star formation. 
The accumulating evi dence in favor of a u niversal initial 
stellar mass function ()Bastian et al.|[2010h indicates that 
star formation on small scales, i.e., within star-forming 
clumps and cores within molecular clouds, is largely de- 
coupled from the ISM properties on larger scales. In 
particular, observations show that the average star for- 
mation efficiency per free-fall time is ~ 0.005 — 0.01, inde- 
pendent of scale, once the densities of mo lecular clouds 
are reached ( e.g.. iKrumholz fc Tanl 120071 but see also 
IMurravl[201ltlFeldmann fc Gnedinll20ilL A natural ap- 
proach is therefore to couple the star formation on scales 
of individual molecular clouds directly to the density of 
molecular hydrogen, assuming that the formation and 
destruction of H 2 can be modeled reliably. 

Consequently, the approach we use in our simulations 
is to relate the SFR to the H 2 den sity on small scales (~ 
60 pc) via the following equation (|Gnedin et al.l l2009): 

P* =csfr h 2 e aX * a , (1) 

TSFR 

where p* is the instantaneous SFR density, pu is the 
hydrogen mass density, /h 2 is the H2 fraction, esfr 
and t sfr denote the star formation efficiency and the 
star formation time-scale (see equation [5] below), respec- 
tively. The exponential factor models intrinsic scatter 
with width a (we discuss intrinsic scatter in £13-21 oth- 
erwise we assume a = 0). A is a Gaussian random 

1 2 

variable with mean and variance 1. The factor e~^ a 
ensures that the same amount of gas is converted into 
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specifies the strength of the interstellar radiation field in 
units of the radiation field of the Milky W ay at 1000A : 
Jmw = 10 6 p hotons cm" 2 s" 1 ster" 1 eV" 1 (|Drainelll97i 
iMathis et afl ll983). We continued one of our simulations 
(Z/Zq = 1, Umw = 1) f° r additional 400 Myr and found 
no significant changes in the Esfr — Eh 2 relation. This 
indicates that the predictions of our simulations should 
also hold for redshifts z < 3, at least unless/until ISM 
properties change radically. In total, we ran a set of 
18 simulations (including 2 runs for a resolution study) 
in order to explore the effect of varying metallicity, ra- 
diation field and density threshold on the Esfr — Eh 2 
relation (see Table [T]). 

The molecular hydrogen fraction /h 2 is computed self- 
consistently, including a chemical network comprised of 
6 species and radiative tra nsfer of the UV continuum and 
the Lyman- Werner bands (|Gnedin fc Kravts"ovll201lD . If 
the average density in a simulation cell is smaller than 
the density typical for molecular clouds, we have to revert 
to a 'subgrid' interpretation of the H2 fraction. In this 
case, we assume that the fraction /h 2 corresponds to the 
(mass) fraction of hydrogen in individual (unresolved) 
molecular clouds. 

Consequently, the star formation timescale is given 
by the minimum of (1) the free-fall time corresponding 
to the average density in the cell and (2) the free-fall 
time corresponding to the minimum density of molecu- 
lar clouds that form stars n c , i.e., 

tsfr = min[r ff (n H ),r ff (n c )]. (2) 
We stress that for densities smaller than n c the relation 
between SFR and H2 abundance is linear, while it be- 
comes non-linear for larger densities, because rg(riH) oc 

— 1/2 

n H . A non-linear steepening of the Esfr — Eh 2 rela- 
tion at Sh 2 > 100M(7) pc~ 2 is motiv ated by theoretical 
studies (e.g., iKrumholz et al.l l2009 h). but not yet con- 
firmed by observations. We therefore explore the case in 
which n c = 50 cm -3 , i.e., close to the typical average 
density of molecular clouds (~ 100 — 200 cm -3 ), but also 
discuss the possibility of much larger thresholds such as 
n c = 10 3 cm~ 3 and n c — 10 6 cm' 3 . Since our simula- 
tions do not capture densities of > 10 5 cm -3 , a threshold 
above this value effectively corresponds to a fully linear 
SFR - H2 relation on small scales. 

Instantaneous SFRs are computed directly using equa- 
tions (fTJ) and @. Our simulations use csfr = 0.005. 
This val ue, which is consistent with small-scale obser- 
vations (jKrumholz fc TarJl200"l . ensures that the nor- 
malization of the Esfr — Ejj 2 relation on kpc scales is 
similar in simulations and observations. Time-averaged 
SFRs (over time T) are calculated by counting the num- 
ber of stars in a cube of given scale with ages below T. 
Unless otherwise noted, we use T = 20 Myr, but we have 
explicitly checked that our results do not change signifi- 
cantly if larger averaging times are used (up to T — 200 
Myr). SFR estimates based on observations of UV lumi- 
nosities in the wavelength rang e 1250-2800 A correspond 
to an averaging time of ~ 100 Myr, which based on neb- 
ular emission lines, such as Ha, typically correspond to 
T ~ 10 Myr, and estimates based on the FIR continuum 
(e.g., 24 /im) correspond to a r ange (~ 10 - 100 M yr) of 
averaging time scales; see, e.g.. iKennicuttl (|1998al ). 

The slope and intercept of the Esfr — Ejj 2 relation 
are obtained with a bisector regression in log-log space 



(|Isobe et alj|1990l ). Although the use of bisector r egres- 
sion cannot be rigorou sly justified in general (see iKellvl 
l2007t iHogg et aLll2010ft , the bisector method is sufficient 
for our purposes as we perform regression on tightly cor- 
related data without error bars. We estimate scatter 
about the best-fit relation, as the root mean square of 
loglO of the spatially averaged star formation rate den- 
sity relative to its value on the regression line with the 
same density, see also equation (jBlj) . We estimate errors 
for the slope, inte rcept, and s catter using the standard 
bootstrap method (|Efronlll979|) with a sample size of 200. 

3. RESULTS 

3.1. Dependence on metallicity and UV field 

In Fig. [1] we plot and compare the Esfr — Eh 2 re- 
lation for (1) solar metallicity and Umw — 0.1, and (2) 
Z/Zq = 0.1 and Umw — 100. Measured over the range 
10 < £h 2 /M q pc~ 2 < 100, the slope of the relation in 
case (1) is ~ 1.14 ± 0.02, the SFR at a surface den- 
sity E H2 = 10 M pc~ 2 is 0.016 M Q yr- 1 kpc" 2 and 
the scatter of log 10 SFR around the best fit is 0.10 dex. 
The slope is slightly steeper tha n that derived fro m CO 
measurements (~ 0.96 ± 0.07; iBigiel et al.ll2008t ). As 
anticipated, the choice esfr = 0.005 leads to a nor- 
malization of the simulated Esfr — Eh 2 relation that is 
close to what is found in observations, once observational 
data are mapped to the same initial stellar mass function 
(IMF). In case (2), the slope is significantly steeper ~ 1.4, 
the SFR at a surface density Eh 2 = 10 M Q pc~ 2 higher 
(0.023 M Q yr" 1 kpc" 2 ) and the scatter is larger (0.18 
dex). 

In Fig. [5] we plot the slope, intercept, and scatter of 
the Esfr — Eh 2 relation, spatially averaged over 1 kpc, 
for a grid of environmental parameters. The figure shows 
that the slope, intercept and scatter are systematically 
changing as a function of Z and Umw- A bi-parametric 
regression (using Z and Umw as independent variables) 
captures the change in slope and scatter very well. The 
regression parameters are given in the legend of the fig- 
ure. A word on the terminology. We refer to the scale at 
which equations (JTJ) and ((2]) are applied as "small scales" 
(~ 60 pc in our simulations). By contrast, we refer to the 
scales on which the slope and intercept of the Esfr — Eh 2 
are measured as "large scales" (~ kpc in this study). The 
spatial density on scale I is given by the amount of mass 
within a cube of size I. In order to convert from a spatial 
to a surface density, we multiply the spatial density by /. 
We do not use the surface density on the smallest (~ 60 
pc) scales in order to avoid underestimating it. 

Fig. [2] shows that the use of time-averaged SFRs intro- 
duces the dominant amount of scatter in the Esfr - Eh 2 
relation on large scales. Specifically, as the first two pan- 
els in the rightmost column demonstrate, the scatter in 
the relation is significantly larger (~ 0.1 — 0.2 dex) if 
SFRs are time averaged, compared with the case that 
instantaneous SFRs are used (~ 0.05 — 0.12 dex). Time 
averaging creates scatter because H2 surface densities are 
measured instantaneously while the SFRs are averaged 
over some past time interval. 

However, the use of time-averaged SFRs is not the 
only source of scatter. The important point to realize 
is that equation (JlJ depends both on the H2 density 
(ph 2 = Ph/h 2 ) an d tne hydrogen density p H (via t S fr)- 
Hence, on small scales, a scatter in the hydrogen den- 
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Fig. 1. — Ssfr — Sjj 2 relation on the kpc scale. Left panel: Z = Zq, Umw = 0.1. Right panel: Z = 0.1 Zq, Umw = 100. Stellar and 
H2 masses are measured within cubical cells of I = 1 kpc box length. Surface densities are estimated by dividing each mass by I 2 . SFRs 
are averaged over 20 Myr. The simulation results are shown by the red shaded region (two-dimensional histogram of all simulation cells) 
and by crosses (a random sample of 50 simulation cells with surface density in the range 10 < E H2 /M pc- 2 < 100). The solid green line 
is the result of a bisector regression of all kpc-sized simulation cells with 10 < Sh 2 /Mqpc -2 < 100 and 0.01 < SFR/MQyr - 1 kpc~ 2 < 1. 
The regression parameters, slope and intercept, are shown on the top left. Also shown (in parentheses) are the regression errors, estimated 
via boo tstrapping. The black circles and triangles correspond to the no rmal spiral and sta r bursting sample, respectively, of IKennicutU 
(1998b). The solid black line is the average Ssfr — ^H 2 relation found in Bigicl ct al. (200<|). The Esfr — Sjj 2 relation in the right panel 
has a steeper slope, a higher normalization, and a larger scatter than in the left panel. 



sity at fixed H2 density translates into a scatter of SFR 
at fixed H2 surface density. The value of the threshold 
n c affects this type of scatter in a crucial way. If n c 
is very large (much larger than the peak in the mass- 
weighted distribution function of molecular hydrogen) 
then the SFR does not depend explicitly on pu (since 
t sfr = Tg(n c )) and, consequently, no scatter is gener- 
ated. Similarly, for hydrogen densities above a certain 
limit, let us call it nf m 5 , the gas is fully molecular and, 
hence , pn and pn 2 are 1:1 related (see, e.g.. lGnedin et all 
2009). If riH > rif m , no scatter is produced on the level 
of a single cell, but scatter can still arise on larger scales 
as cells with different properties are added. To clarify 
this point, let us assume that we add the SFRs and H2 
densities from, e.g., two cells A and B. First, let cell A 
have a density below n c and cell B a density above rif m . 
Second, let us redistribute the hydrogen and H2 masses 
such that both cells have a density below n c (this might 
not be possible in all cases). Although in both cases the 
H2 density is the same, the SFRs are higher in the first 
case. 

The mechanism that we have just described explains 
the existence of scatter, provided n c is sufficiently low 
(see the third column of Fig. [2J. However, we have not 
discussed why there is a trend of scatter with Z and 
Umw- The origin of this trend can be understood from 
the bottom panels of Fig. [3l where we show a histogram 
of the small-scale hydrogen density (H2 mass-weighted) 
parametrized by the large-scale H2 surface density. The 
figure shows that the fraction of H2 mass that is in cells 
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with hydrogen densities above a given threshold (in the 
range of ~ 10 — 100 cm -3 ) increases with decreasing Z 
and increasing UV. Hence, more of the H2 mass partic- 
ipates in producing scatter and the overall scatter in- 
creases. 

Fig. [2] also shows that there is a dependence of the 
intercept of the Ssfr — £h 2 relation on Z and Umw, 
provided n c is sufficiently small (see the middle panel in 
the top and middle rows). How do we understand this 
result? As we just pointed out, an increase in the ra- 
diation field and/or a decrease in the metallicity shifts 
the peak of the mass weighted H2 density distribution 
towards higher densities (see the bottom row of Fig. [3]). 
Specifically, the figure shows that for Umw = 0.1 and 
Z/Zq = 1, cells with hydrogen densities in the range 
of ~ 3 — 100 cm -3 contain most of the molecular hy- 
drogen for the considered range of H2 surface densities. 
Hence, only cells in the range of ~ 3 — 100 cm -3 con- 
tribute significantly to the SFR. On the other hand, if 
Umw = 100 and Z/Zq = 0.1, only cells with hydrogen 
densities in the range of 100 — 500 cm -3 contribute to 
star formation. Hence, a low-metallicity, high UV disk 
will only keep "pockets" of H2 in high density regions, 
while in a high-metallicity, low UV disk H2 is present 
even in much lower density regions. Consequently, for a 
given large-scale H2 surface density, more of the H2 sits 
at high densities in a low Z, high UV field, galaxy com- 
pared to a high Z, low UV field galaxy. Furthermore, 
in the regime in which riu > n c and the hydrogen gas is 
(close to) fully molecular the SFRs scale as oc n^ 5 . This 
non linearity then amounts to a higher SFR (and hence 
intercept) at fixed large-scale H2 surface density in a low 
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Fig. 2. — Dependence of slope, intercept, and scatter on metallicity and interstellar radiation field. The first row shows (from 
left to right) the slope, the intercept, and the scatter of the Esfr — Eh 2 relation (averaged over kpc scales) as a function of the 
radiation field, Umw (a;-axis), and for different metallicities (Z/Z® = 1 (black squares), 0.3 (red downward-pointing triangle), 
and 0.1 (blue upward-pointing triangle)). SFRs are averaged over 20 Myr and the minimum cloud density is n c = 50 cm -3 . 
Slope, intercept and scatter are computed from a bisector regression of the Esfr — Eh 2 relation as described in the caption of 
Fig- HI To highlight the trends with Z and Umw, we also performed a two-parametric regression of slope, intercept, and scatter 
as a function of Z and Umw (regression equation and parameters and the square of the correlation coefficient are shown at the 
top of each panel; xi = log 10 Umw, %2 = log 10 Z). The black, red, and blue solid lines (from bottom to top) show the results 
of the bi-parametric regression for the choices Z/Zq = 1, 0.3, and 0.1, respectively. The middle row shows the same quantities 
as the top row, but for instantaneous SFRs. The bottom row shows again the same quantities, but for a larger threshold 
density n c . Specifically, the filled symbols and lines use instantaneous SFRs and n c = 1000 cm" 3 , while the empty symbols use 
time-averaged SFRs and n c — 1000 cm -3 (small magenta symbols) and n c = 10 6 cm -3 (large cyan symbols), respectively. We 
note that whenever n c > 50 cm -3 the star formation efficiencies are reduced by -\/n c /50 (see equations {TJ and ©) in order to 
ensure the correct normalization of the Esfr — Eh 2 relation. 
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Fig. 3. — Distribution of neutral hydrogen mass (top row) and 
H2 mass (bottom row) as a function the small-scale (~ 60 pc) 
hydrogen density tijj (ordinate) for several large-scale (1 kpc) H2 
surface densities (abscissa). Each vertical (H2 surface density) bin 
is normalized to the total mass of neutral hydrogen (top row) or 
H2 (bottom row) in the bin. The grade of shading indicates the 
mass fraction (on a linear scale). A black pixel contains 100% of 
the neutral hydrogen or H2, respectively, in the given H2 surface 
density bin; white corresponds to 0%. The left column corresponds 
to a galaxy with high metallicity and low radiation field (Z = 
Zq, Umw = 0.1), while the right column is for a low- metallicity, 
high radiation field galaxy (Z = 0.1 Zq, Umw = 100). Mean and 
standard deviation of the density distributions of neutral hydrogen 
(top row) or H2 (bottom row) are indicated by green crosses and 
error bars. The thin, red horizontal line indicates the threshold 
n c = 50 cm -3 . The thick, black horizontal line near the bottom of 
each panel shows the range over which slope, intercept, and scatter 
of the Esfr — Sjj 2 relation are computed, see Fig. [2] 

Z, high UV field, galaxy compared to a high Z, low UV 
field galaxy. In other words, the large-scale SFRs depend 
not only on the large-scale H2 surface densities, but also 
on the distribution function of jih on small scales. 

A related mechanism leads to a dependence of the slope 
on Z and Umw (left panel in the top and middle rows of 
FigJ2]). The panels in the top row of Fig. [3] show that 
the typical hydrogen densities of cells that contribute 
to a given H2 surface density £h 2 increase with £h 2 - 
This demonstrates that the density structure of disks at 
a given H2 surface density changes with metallicity and 
radiation field of the ISM. The bottom panels in Fig. [3] 
show that this trend remains (although somewhat weak- 
ened) if the hydrogen density distribution is weighted by 
H2 mass. It also shows, that the effect of an increasing 
hydrogen density with large-scale H2 surface density is 
stronger for a high Z, low Umw galaxy than for a low 
Z, high Pmw galaxy. However, in the former case, most 
of the H2 is in cells with densities nn < n c and conse- 
quently the SFR density is still proportional to the total 
H2 density, i.e., slope 1. In the latter case, however, 
this increase is important. Let us see why. From (fTJ) 
and ([2]) it is clear that the surface density of the SFR 



is proportional to J2i n H 



where the sum is over 



all cells within the given line-of-sight cylinder and a is 
either (if nn < n c ) or 0.5 (if nn > n c ). The H2 sur- 
face density is proportional to X)i n H 2 ,i- An increase in 
the H2 surface density Sh 2 — > 7^h 2 can be achieved in 
several ways. If the length of the cylinder increases, then 
the surface density of the SFR increases proportional to 
Eh 2 - If, however, the density structure changes (in the 
simplest case via nn ~> 7^h), then the surface density 
of the SFR increases by j 1+a (assuming the gas is fully 
molecular). Obviously, if a = (as for nn < n c ) the 
predicted large-scale slope of the Esfr — £h 2 relation is, 
as expected, linear. However, it lies between 1 and 1 + a 
if there is a mix of cells with densities below and above 
n c . In addition, if the density distribution changes in a 
more complicated manner with Sh 2 , it is also possible 
to obtain "rolling" slopes or even large-scale slopes that 
are steeper than 1.5. We conclude that the slope of the 
Ssfr — Sh 2 relation on kpc scales can vary and depends 
on the H2 density distribution as a function of th e large- 
scale H2 surface density, see also iKravtsovl ()2003f ) . 

While the time averaging of the SFRs generates most 
of the scatter in the Ssfr — ^h 2 relation, the trends of 
slope and scatter with Z and Umw are largely driven 
by the non-linear coupling between SFR and H2 density. 
This can be clearly seen in the last row of Fig. [2] If n c — 
10 6 cm -3 , the slope of the Ssfr ^ £h 2 relation changes 
only between 1.03 (Z/Zq — 1, Umw = 0.1) and 1.16 
(Z/Zq = 0.1, Umw = 100) and the scatter increases only 
from 0.09 dex to 0.12 dex. We discuss the dependence of 
the scatter on ISM properties further in the next section. 

If SFRs are measured instantaneously and the small- 
scale relation between star formation rate density and H2 
density is linear (i.e., n c is large), then the slope reduces 
to exactly unity, and any dependence of the intercept 
on metallicity or radiation field is eliminated and the 
scatter vanishes (at least as long as there are no other 
sources of scatter, see Ej3.2[) . The reason lies in the fact 
that spatial averaging (which is a linear operation) over 
a linear relation between SFR and H2 density on small 
scales results again in a linear Ssfr — Sh 2 relation on 
larger scales. 

We conclude that slope, intercept, and scatter of 
the Esfr — Sh 2 relation averaged on kpc scales can 
change systematically with metallicity and radiation 
field. While our quantitative predictions likely depend 
on the assumed star formation model 6 and the model 
for H2 formation and shielding the existence of a "non- 
universality" of the large-scale Ssfr — £h 2 relation is 
a rather generic outcome whenever there is a non-linear 
relation between SFR and H 2 density on small scales. 

3.2. Dependence on averaging scale 

Observational studies show that the Es fr ~ ^Hi re- 
lation has larger scatt er on sm aller scales (|Verlev et aLl 
120101 lOnodera et ail l2010t iDanielson et all 120111) . 
Specifically, recent observations of CO, Ha, and 24^m 
emission in M33 have been used to argue that the 
^SFR — Sh 2 relation "breaks down" on a scale of ~ 100 
pc. It has been suggested that the drifting of newly 
formed star clusters or the difference in evolutionary 
stages of molecular clouds / s tar clusters could be re- 
sponsible (jOnodera et al.ll2010D . Given the limited range 

6 For instance, as we show in Fig. [2] they vary with n c . 
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Fig. 4. — Dependence of the scatter in the Ssfr — Sh 2 re- 
lation on the averaging scale. The scatter in the SgpR — £h 2 
relation has been derived assuming either (1) time-averaged SFRs 
(20 Myr) and a minimum cloud density of n c = 50 cm -3 (black 
squares), (2) time-averaged SFRs and n c = 10 6 cm -3 (red tri- 
angles), or (3) instantaneous SFRs and n c = 50 cm -3 (blue dia- 
monds). The scatter is computed with an ordinary least squares 
regression of the SFR as a function of H2 surface density in the 
range 10 < Sh 2 /Mqpc — 2 < 100. Errors are computed via 
bootstrapping. Empty symbols correspond to solar metallicity 
and Umw = 0.1, while filled symbols refer to Z = 0.1 Zq and 
Umw = 100. Simulations with intermediate values of metallicity 
and UV field lie in between. 

of measured gas surface densities, it is plausible that 
this "break-down" is merely a manifestation of a very 
large scatter that ma y arise from a variety of sources. In 
lOnodera et al.l (|2010t ). the studied range of surface den- 
sities is approximately lM Q pc -2 < £h 2 < 10Af©pc~ 2 . 
By contrast, the average gas surface density of GMC 
measured on s cales of a few tens of p arsec in M33 is 
~ 120M Q pc -2 (|Rosolowskv et alj|200ai, roughly similar 
to that found in the Milky Way (|Solomon et al.l 119871: 
iHever et al.ll2009( ). Hence, H 2 gas surface densities be- 
low < 10M Q /pc 2 measured on 100 pc scales must cor- 
respond to the outskirts of GMCs, not to GMCs them- 
selves. This by itself may be responsible for a substan- 
tial fraction of the measured scatter. We, instead, will 
focus on a range of 10 times larger H 2 surface densities, 
namely 10M Q /pc 2 < E H2 < 100M o /pc 2 for scales from 
kpc down to 100 pc. Instead of aiming at a precise quan- 
titative comparison with observations (that are not yet 
available), we discuss qualitative predictions of our simu- 
lations, their limitations, and the resulting implications. 

In Fig. [4] we show the scatter of the Ssfr — £h 2 re- 
lation as a function of scale from 1 kpc down to ~ 100 
pc. The first thing we notice is that the scatter due to 
time-averaging alone (red lines and triangles) increases 
with decreasing scale, while the scatter solely due to the 
threshold density n c (blue lines and diamonds) remains 
roughly scale-independent. The origins of the different 
types of scatter have been discussed in the last section. 



The figure shows that the scatter on all scales is primarily 
caused by the time averaging of the SFR. This is not an 
artifact of the particular SFR averaging time used. We 
varied the SFR averaging time scales between 20 and 
200 Myr and found no substantial change in the amount 
of scatter, as long as low SFR outliers (> 3 sigma) are 
excluded 7 . 

A straightforward comparison between our predic- 
tions and observations |Verlev et al.ll2010HOnodera et al.1 
I2010t ISchruba et al.ll2010D is difficult as none of the pub- 
lished observations explicitly quantify the scatter as a 
function of scale. However, the scatter predicted from 
our simulations (~ 0.4 — 0.6 dex on ~ 100 pc, ~ 0.1 — 0.3 
on ~ kpc scales) seems to be of similar order as, e.g., the 
observed scatter in Fig. 4 of IVerlev et al.l (|2010f) and the 
scatter in Fig. 3 of ISchruba et al.1 (|2010f ). One should 
keep in mind though, that the scatter in the literature 
may include contributions from obse rvational uncertain - 
ties. Also, the standard deviation in lVerlev et al.l (2010) 
is computed on iteratively sigma-clipped data, which 
likely underestimates the scatter. Even more problem- 
atic is that all mentioned observational studies derive H 2 
masses from CO measurements, which may introduce ad- 
ditional scatter. A detailed modeling of the H2-CO con- 
version factor which would be required for a fair compar- 
ison is beyond the scope of this paper. 

Fig. |4] also shows that the scatter due to the SFR time 
averaging depends to a small extent on Uuw and Z. 
This result can be rephrased in terms of a duty frac- 
tion, which we define as the fraction of time (the rele- 
vant time scale is the SFR averaging time) during which 
the H 2 density within the considered cell is close to its 
time-averaged value. A duty fraction of unity does not 
introduce scatter on small scales, as it means that within 
the SFR averaging time the H 2 content within the cell 
remains constant. Fig. [4] then shows that the duty frac- 
tion decreases with increasing radiation field/decreasing 
metallicity, hence leading to larger scatter. One interpre- 
tation of the reduced duty fraction is that stronger Umw 
and/or lower Z reduce the life times of molecular clouds. 
An alternative possibility is that molecular clouds live as 
long as before, but molecular cloud formation is rarer. 

Another potential contributor to the scatter on small 
scales is the velocity spread of young stellar clusters and 
the stars within the cluster. This effect is not modeled 
adequately in the simulation because we do not resolve 
individual cluster members, but rather obtain one 'star 
particle' for each cluster that initially moves with the 
average velocity of the gas. On scales > 100 pc this 
effect plays only a small role presumably, as the typical 
distance that stars travel within 20 Myr is of the order 
of ~ 100 pc (assuming an rms velocity of ~ 5 km s _1 ). 
Some scatter on large scales may arise from high- velocity 
run-away stars (|Blaauwlll96"lt IStondll99lh . 

The decline of the scatter with increasing averaging 
scale is obviously related to the spatial averaging over a 
larger number of resolution elements N rcs . Naively, we 
would expect a scaling proportional to l/^/N les , where 
Arcs oc Z 3 if the H 2 is filling the volume relatively uni- 
formly, or oc I 2 if the H 2 is confined to a disk. However, 
the scale dependence shown in Fig. [4] seems to be much 

7 If not excluded, these outliers do increase the scatter somewhat 
(by ~ 0.1 dex) when averaging over 200 Myr instead of 20 Myr. 
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Fig. 5. — Propagation of scatter in the Sspr — Eh 2 relation from 
small to large scales. A scatter in log 10 (SFR[Moyr~ 1 ]) of (from 
bottom to top) 0.1 (blue lines), 0.25 (cyan lines), 0.5 (magenta 
lines), and 1 (red lines) is inserted at the resolution scale (65 pc) of 
the simulation. The black line at the bottom shows a power law: 
scatter <x scale -0,5 . The scatter in the EsFR — Eu 2 relation is 
measured in the range 10 < Sh 2 /Mqpc — 2 < 100. Instantaneous 
SFRs and a minimum cloud density of 10 6 cm -3 are used in order 
to suppress other sources of scatter. Filled (empty) symbols corre- 
spond to the simulation with Z = 0.1 Zq, Umw = 100 (Z = Zq, 
Umw = 0.1). 



shallower. 

In order to understand this better, we now include a 
log-normal (intrinsic) scatter of the SFRs 8 in equation 
J]}, with cr/lnlO ranging from 0.1 to 1 dcx on the 65 
pc scale. Fig. [5] shows that this scatter decreases with 
increasing averaging scale roughly as a power law 07 oc 
l~ a with exponent a ~ 0.5. 

As we discuss in more detail in Appendix [Bl the rather 
gradual decline is caused by both the finite width a p 
of the H2 density distribution and the geometrical ar- 
rangement of the molecular hydroge n. In particular, we 
find that a sensible value a p ~ 1.5 (M cKee fc Ostrikerl 
|2007|) and a two-dimensional (2D, disk) arrangement of 
the molecular hydrogen naturally leads to a ~ 0.5. On 
the other hand, a purely one-dimensional (ID) arrange- 
ment of the H2 is only consistent with a ~ 0.5 if the 
width in the density distribution is very small (<C 1), 
while a three-dimensional (3D) configuration would re- 
quire a p w 2.5. We have also tested that the exponent 
decreases if the inserted intrinsic scatter becomes large, 
i.e., u >■ 1. In fact, the top curves (a — 2.3) in Fig. [5] 
clearly show this flattening. Since the H2 density dis- 
tribution depends on the metallicity and radiation field 
(see Fig. [2]), a should depend on it, too. With the data 

8 In the remainder of this section we refer to instantaneous SFRs. 



points shown in Fig. [5l we obtain a = 0.52 ± 0.04 for 
Z = 0.1, Umw = 100, and a = 0.43 ± 0.04 for Z = 1, 
Umw = 0.1. 

Let us now consider the case in which a scatter 07 is 
inserted on scale I (we assume a set of discrete scales that 
change by a factor 2). If the different scatter contribu- 
tions add in quadrature, the total scatter 07 on a given 
scale I is simply given by 



of- 



~2 



1 *2 



1 

42S 



^ 2 



(3) 
(4) 



With knowledge of a, this equation allows the compu- 
tation of the amount of scatter 07 that is introduced on 
scale I from the measurement of the scatter on scales 
I and 1/2. Presumably, different physical mechanisms 
may introduce different amounts of scatter on different 
scales. Studying the scale dependence of the scatter may 
therefore be helpful to uncover the responsible physical 
mechanism(s). 

4. DISCUSSION 

The scatter in the Esfr — £h 2 relation on the scale of 
~ 100 pc has been attributed to the evol ution of molec- 
ular clouds over their life time (see, e.g.. lOnodera et al.l 
l2010HSchruba et aJJ[20Toh . In this picture, young molec- 
ular clouds have not yet formed stars, but contain 
large amounts of H2 and hence fall "below" the aver- 
age SgpR — Eh 2 relation. On the other hand, clouds that 
are near the end of their lives are heavily star forming 
and/or have lost some fraction of their molecular hy- 
drogen, hence they lie "above" the relation. This picture 
cannot be reconciled with an H2-based star formation law 
of the form of equation (JXJ) as long as the gas consump- 
tion time scale tsfr/esfr is treated as a constant. Hence, 
this explanation of the scatter in the Esfr — ^h 2 rela- 
tion implies that csfr/tsfr has to be a time-dependent 
quantity. If tsfr is approximately constant, then the 
star formation efficiency will need to change ove r the life 
time of a molecular cloud (e.g., iMurravl 120111 but sec 
iFeldmann fc Gnedidl2l)ll . 

Our interpretation is different. We show that a large 
amount of scatter in the Esfr — ^h 2 relation can be ex- 
plained by the fact that observations do not measure the 
instantaneous rate of star formation, but rather the num- 
ber of stars that formed within a finite time interval in 
the past. Our numerical models predict that the scatter 
seen on scales of ~ 100 pc should be small (< 0.1 dex) if 
SFRs are measured instantaneously. By contrast, if the 
scatter in the relation is mainly due to an evolving star 
formation efficiency, the scatter on 100 pc scales should 
not be strongly diminished if (close to) instantaneous 
SFRs are used. 

We note that the particular small-scale model of star 
formation used in our simulations is based directly on 
the H2 density (equation [I]) . Although this mode l has 
been motivated analytically (Krumholz et al. 2009 b]) and 



is wi dely used in numerical simulations (Gned in et al 



2009) or, more recently, semi-analytic models (|Fu et al 
2010T ) . it is important to verify how accurately it de- 



scribes reality. A potential shortcoming of equation ([T]) 
is that it assumes that a fixed mass fraction of molecu- 
lar hydrogen is available/eligible for star formation. It 
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has been known for a while now that star formation in 
molecular clouds occurs preferentially in t he region of 
high gas density (n ^ 10 4 cm" 3 , see e.g., lLadal[l99l: 
iGao fc Solomon|[200llLada et al.ll2010l) . Hence, the SFR 
should be strongly dependent on the density probability 
distribution function (pdf ) of the gas and not necessarily 
on the total mass of molecular hydrogen alone. 

The density pdf can be expected to depend on details 
of gas t hermodynamics or poten ti al feedback mechanisms 
(e.g., IWada fc Norman! hOQlt iRobertson fc Kravtsovl 
2008), even in a picture in which the turbulence in the 
ISM is ma inly driven by large-scale gravitational mo - 
tions (e.g., Wada fc Normanl[200l iTasker fc Tanl 12009). 
Hence, if the SFR is in fact regulated by the amount of 
high-density gas (and not by H2), we can expect to see 
differences in the SFR on ~ kpc patches as a function of 
Z and Umw, even for the same star formation prescrip- 
tion on small scales. 

In addition, the total amount of H2 may depend on Z 
and UMW. For instance, the total amount of H2 in the 
simulated volume changes by a factor of ~ 2 — 3 when 
metallicities and radiation fields are varied in the range 
Z/Zq =0.1-1 and Uuw = 0.1 - 100. This change in the 
H2 mass alone should induce a galaxy-by-galaxy scatter 
on the order of ~ 0.2 dex. A study of the scatter of the 
Esfr — Sh 2 relation as a function of scale and in regions 
with different metallicities and UV radiation fields may 
thus give us, at least in principle, a means to test this 
picture. 

5. CONCLUSIONS 
5.1. The scatter in the Esfr — Eh 2 relation 

Our simulations identify and quantify two important 
sources of scatter. The first type of scatter is only present 
if the small-scale star formation relation is non linear and 
arises due to the scatter in the H 2 density at fixed gas 
density. This scatter is relatively independent of spatial 
scale and amounts to (at most) ~ 0.1 dex. The second 
type of scatter is due to the fact that observations of H2 
(or CO), which measure the instantaneous gas surface 
density, are combined with measurements of SFRs that 
are averaged over the past tens of megayears. This scat- 
ter does not strongly depend on the averaging timescale 
provided it is longer than both the H2 formation and de- 
struction time, i.e., it exceeds ~ 10 Myr 9 . For shorter av- 
eraging times, the scatter should decrease and will even- 
tually be dominated by other scatter sources. This type 
of scatter is strongly dependent on spatial scale. It varies 
between ~ 0.1 — 0.2 dex (in log 10 SFR) on kpc scales and 
~ 0.4 — 0.6 dex on ~ 100 pc scales. 

Our numerical experiments predict that intrinsic scat- 
ter scales with averaging scale I approximately as cx Z~ - 5 . 
This relatively shallow scaling is primarily caused by the 
finite width of the H2 density distribution and the ar- 
rangement of the molecular hydrogen in a 2-D (disk) 
configuration. 

We note that our simulations provide only a lower limit 
on the expected scatter as a function of scale, because 
some sources of scatter (e.g., the velocity dispersion of 
star clusters and their member stars, or scatter in the 

9 Assuming typical conditions rajj ~ 50 cm -3 , T ~ 80 K of 
the cold neutral medium and tn 2 = nn/hn 2 = l/(Rdnn) with 
R d = 6 X 10" 18 T 1 / 2 cm 3 s" 1 jDraine fc Bertoldill996h . 



CO to H2 conversion) are not included in our numeri- 
cal modeling. A precise observational determination of 
the scatter-scale relation, possibly even as a function of 
ISM environment, and the comparison with theoretical 
predictions, such as the one presented in this paper, may 
thus help to identify the physical processes responsible 
for creating the scatter. Consequently, we argue that 
the scale dependence of the scatter in the Esfr — Eh 2 
relation could become an important diagnostic tool in 
determining the underlying connection between star for- 
mation and H2 density. 

5.2. The environmental variation of the Esfr — Ejj 2 

relation 

We have shown that even if the SFR is tightly cou- 
pled to the H 2 density on small scales (see equation [1]) , 
the Esfr — Ejj 2 relation can vary systematically with 
metallicity and interstellar radiation field in the studied 
surface density range 10 < Eh 2 / 'M pc~ 2 < 100, when 
averaged on ~ kpc scales. 

In particular, the super-linear slope of the Esfr — Eh 2 
relation depends on the actual H2 density distribution 
and on the existence of a non-linear scaling between SFR 
and H2 density. The underlying reason for a slope steeper 
than unity is that the peak of the H2 density distribution 
changes systematically with large-scale surface density. 
At larger ~ kpc averaged surface densities, more of the 
molecular gas sits at higher densities, which, due to the 
non-linear scaling between SFR and density, leads to the 
super linear steepening of the Esfr — Eh 2 relation. 

Similarly, the systematic change in the slope with 
metallicity of the ISM and the interstellar radiation field 
is a reflection of the change in the H2 density distribution. 
For example, in a low metallicity and/or strong radiation 
field environment, the HI to H2 transition takes places at 
significantly higher densities and, consequently, a larger 
fraction of the H2 mass contributes super-linearly to the 
SFR. In addition, this implies more star formation at a 
given H2 surface density and hence changes the intercept 
of the Esfr — Eh 2 relation. 

The scatter in the Esfr — Eh 2 relation also shows a 
systematic trend with Z and Umw- The precise value 
of the scatter and the amount it changes with Z and 
Umw depends on (1) the assumed density threshold, n Cl 
above which the SFR scales super-linearly with density, 
and (2) amount of time over which observed SFRs are 
time averaged. The scatter varies between ~ 0.05 dex 
{Z/Zq = 1, Umw = 0.1) and - 0.12 dex (Z/Z Q = 0.1, 
Umw — 100) if n c = 50 cm -3 and SFRs are measured 
instantaneously. The trend with Z and Umw is mainly 
caused by changes in the H2 and total hydrogen density 
distributions. On the other hand, if n c is large (> 10 4 
cm~ 3 ) and the scatter is generated by the time averaging, 
then the scatter changes only weakly with Z and Umw 
(~ 0.09 dex versus ~ 0.13 dex). 

We note that in order to observe a significant change in 
the Esfr - Eh 2 relation, metallicities < 0.3 Zq and inter- 
stellar UV fields Umw > 10 are required. Star-forming 
galaxies at high redshifts should therefore be the natural 
candidates to test our predictions. 
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APPENDIX 
RESOLUTION STUDY 

In order to test for possible resolution effects we have rerun one of our simulations (Z = 1 a nd Umw = 1) at two 
times better (32 pc) and also two times worse (125 pc) spatial resolution. We show in Fig. I All the scale dependence 
of the scatter in the £sfr — ^h 2 relation for the three different resolutions (32, 65, 125 pc). As in Fig. [U we present 
both the scatter (1) due to SFR time averaging (using a high n c ) and (2) due to the non-linear scaling between SFRs 
and H2 density (using instantaneous SFR). We find that the amount of scatter as a function of scale is similar in each 
of the three simulations and that there is no apparent significant systematic trend with resolution. 

DEPENDENCE OF THE SCATTER ON SPATIAL AVERAGING SCALE 

As discussed in section |3~21 the scatter in the Esfr — ^h 2 relation decreases rather slowly (oc l~ - 5 ) with averaging 
scale I, much slower than expected from a U~ ' 5 scaling, unless the gas is arranged in a ID configuration. The scatter 
(j; in the Ssfr — Sh 2 relation can be formally written as 

= ((logiopT- (log 10 pT)) 2 ), (Bl) 

where p* = y J dxp\(x) is the spatial average of the star formation rate density and the brackets (. . .) denote the 
ensemble average at fixed spatially averaged H2 density p = h J dxp(x). 

The star formation rate density p» is a random field and its value for a given point in space is defined by equation 
©, i.e., 

Pxocpe' 7 *, (B2) 

where p = pn 2 and we assume that p» is proportional to the H2 density, i.e., tsfr is a constant. Here A is a Gaussian 
random field with zero mean and unit variance. 
From equations (|B1[) and (|B2[) . it follows that 07 = aj In 10 on the scatter insertion scale Z m i n , i.e., without spatial 
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Fig. Bl. — Scatter ct; of the SgpR — Sh 2 relation as a function of the averaging scale I for different widths of the log-normal density 
distribution a p and different geometrical arrangements of the molecular hydrogen. The dot-dashed and solid lines show the result of a 
simple numerical experiment (see text), in which a spatially uncorrelated scatter with cr; min = 0.5 is inserted on the scale of ( m i n = 100 pc 
and consequently spatially averaged on larger and larger scales. The dot-dashed lines correspond to a p = 0.1 in three (blue, bottom) and 
one (magenta, top) dimensions, respectively. The red solid lines shows the result for a 2D gas configuration with a p = 0.1 (lower line) and 
Up = 1.5 (upper line). Dashed lines (black) indicate the scaling ex l~ 0,5 J oc l —1 , and tx l~ 15 (from top to bottom) and are shown to guide 
the eye. 

averaging. The reason being that, in this case, — /?*, p = p, (log 10 p) = log 10 p and thus 

= ( ( l0g - M + hr!) - ^M) - (^)) ) (B3) 
= ((log w p+^--{log w p)) \ (B4) 



' aX 
In 10 




(B5) 



Equation (|B1[) can be solved if all cells have the same density p and a <^ 1. In this case epe aX = ep e aX , 
log 10 e aX ~ aX, and 

of = <(log 10 (i^) - (log 10 (^))) 2 ) (B6) 



In 10 



2 



-2 



(X ). (B7) 



2 

If X is spatially uncorrelated then (X(x)X(y)) — S(x — y), (X ) = 1/V, and <ti oc V~°^. 

However, the scaling of a\ may be different if X is spatially correlated, a 3> 1, or if the density pdf is not sharply 
peaked. We can test whether the latter, i.e., the finite width of a log- normal density distribution, provides a quantitative 
explanation for the scaling of ai. To this end, we insert spatially uncorrelated scatter with <r/ min = 0.5 on the scale 
of 100 pc and compute ai as a function of averaging scale for log-normal density distributions of various widths a p 
and for a ID, 2D and 3D configuration of H2. Specifically, for each scale I we first generate (l/l m in) d independent 
density values (d is the assumed dimensionality of the gas configuration) drawn from a log-normal distribution with 
width <7 P and the same number of corresponding star formation rate density values according to equation (|B2j) . We 
then compute the average H 2 density and SFR over the (l/l m i n ) d "cells". This pair of spatially averaged H 2 density 
and star formation rate density constitutes a data point in the Esfr — Sh 2 relation and the scatter of the relation is 
computed from 1000 data points generated in this way. 

Fig. IB1I shows that if a p is sufficiently small (<§; 1), 07 scales as V~ ' 5 as anticipated. Supersonic turbulence 
simulations, however, predict a p = ln(l + fi 2 M 2 ) 1/2 with M ~ 5 - 10 and /3 ~ 0.25 - 0.5 ()McKee fc Ostriker|[2007l) . 
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Thus, dp Ri 1.5 is a more reasonable assumption. In this case, 07 scales approximately oc 
a disk. 



-0.5 



if the gas is arranged in 
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